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Abstract:  The  optimal  sensor  placement  problem,  as  considered  here,  is 
to  select  the  types  and  locations  of  sensors  providing  coverage  at  high- 
value  terrain  locations  while  minimizing  a  specified  cost  function.  The  cost 
function  can  reflect  various  disincentives,  such  as  the  actual  cost  of  the 
sensors,  the  total  number  of  sensors,  and  the  probability  that  the  sensor 
will  be  found  and  disabled  by  hostile  actors.  The  probability  of  detection 
(at  a  certain  probability  of  false  alarm)  is  assumed  to  depend  on  terrain 
conditions  and  obstructions,  and  may  be  arbitrarily  complex.  Two  strate¬ 
gies  are  described  for  finding  the  minimal  number  of  sensors,  and  their 
locations  that  will  satisfy  given  coverage  preferences.  The  first  is  heuristic 
in  nature  and  based  on  placing  sensors  one-by-one  where  the  probability 
of  detection  is  minimal.  This  strategy  offers  a  rapid,  but  suboptimal  solu¬ 
tion.  The  second  strategy  is  based  on  solution  of  the  binary  linear  pro¬ 
gramming  problem.  For  the  case  of  fine  spatial  resolution  that  leads  to 
large  matrix  dimensions,  a  fast  algorithm  for  approximate  solution  of  this 
problem  is  developed.  The  key  features  of  this  study  are:  l)  the  probabilis¬ 
tic  framework  of  sensor  performance,  2)  incorporation  of  the  coverage 
preferences  in  the  placement  strategy,  3)  realistic  modeling  and  incorpora¬ 
tion  of  the  sensors’  probability  of  detection,  4)  multimodal  sensor  support, 
5)  a  strict  formulation  of  the  optimal  coverage  problem,  6)  development  of 
a  fast  algorithm  for  approximate  solution  of  the  binary  linear  program¬ 
ming  problem,  and  7)  introduction  of  a  safe-mode  concept. 
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1  Introduction 

The  problem  of  optimal  sensor  placement  can  be  viewed  as  a  generaliza¬ 
tion  of  the  art  gallery  problem  originally  formulated  as  finding  the  mini¬ 
mal  number  and  locations  of  guards  so  that  each  point  of  a  gallery  would 
be  seen  by  at  least  one  guard.  This  problem  and  its  variations  are  well 
studied  in  computational  geometry,  and  elegant  solution  algorithms  have 
been  developed  (O'Rourke  1987;  Shermer  1992).  Because  the  problem  is 
obviously  identical  to  the  problem  of  camera  placement,  it  is  not  surpris¬ 
ing  that  these  algorithms  are  being  used  for  security  and  surveillance.  The 
problem  can  be  generalized  to  include  constraints  on  the  camera's  range, 
resolution,  and  field  of  vision  (Erdem  and  Sclaroff  2004).  Similar  prob¬ 
lems  are  encountered  in  many  other  areas,  such  as  broadcasting  or  cellular 
tower  placement  (Eidenbenz  2002),  robot  motion  planning  (Elnagar  and 
Lulu  2005;  Gonzalez-Banos  et  al.  1999),  graphics  (McKenna  1987),  and 
computer  vision  and  pattern  recognition  (Mittal  and  Davis  2008;  Wixson 
1994). 

In  this  technical  report,  the  original  art  gallery  problem  is  generalized  to 
include  sensors  of  different  modality  (e.g.,  cameras  and/or  acoustic,  infra¬ 
red,  seismic,  and  magnetic  sensors).  To  describe  and  operate  with  multi¬ 
modal  sensors  systematically,  the  probability  of  detection  (at  a  certain 
probability  of  false  alarm)  is  chosen  as  a  sensor's  performance  measure 
(Wilson  et  al.  2008);  it  is  universal1,  it  directly  relates  to  the  receiver  op¬ 
erating  characteristic  (ROC),  which  is  a  standard  characteristic  of  a  sensor, 
and  it  provides  a  natural  way  to  formulate  coverage  preferences  by  directly 
assigning  the  desired  probability  of  detection  for  each  spatial  point.  (In 
this  technical  report,  a  standard  terminology  of  signal  detection  is  used: 
the  probability  of  detection  Pd  is  the  probability  that  a  sensor  would  detect 

a  source  when  it  is  really  there;  the  probability  of  false  alarm  Pfa  is  the 

probability  that  a  sensor  would  detect  a  source  when  it  is  really  not  there; 
the  probability  of  misdetection  Pmd  =  1  -  Pd  is  the  probability  that  a  sensor 

would  not  detect  a  source  when  it  is  really  there.)  However,  this  brings  the 


1 A  probability  of  detection  can  be  applied  to  any  sensor.  For  example,  an  ideal  camera  has  a  probability 
of  detection  equal  to  one  within  its  field  of  view  and  zero  outside. 


ERDC/CRREL  TR-08-24 


2 


problem  into  probabilistic  space,  which  generates  two  new  aspects  in  com¬ 
parison  with  the  ordinary  art  gallery  problem. 

First,  the  ordinary  art  gallery  problem  and  its  variations  require  that  any 
point  within  a  specified  area  be  seen  by  at  least  one  guard  (or  camera). 

This  requirement  is  unnecessary  in  the  probabilistic  framework  of  sensor 
performance.  Indeed,  to  have  sufficient  coverage  at  a  certain  point  (that  is, 
a  sufficiently  high  probability  of  detection,  for  example,  over  0.95 ),  it  is 
not  necessary  to  require  that  at  least  one  sensor  see  it  (that  is,  that  at  least 
one  sensor  have  a  probability  of  detection  over  0.95  at  this  point).  For  ex¬ 
ample,  two  independent  sensors  can  have  a  probability  of  detection  at  a 
certain  point  equal  to  0.8  to  see  it  with  a  probability  of  0.96 . 

Second,  the  effects  of  signal  propagation  from  the  point  where  a  signal  was 
generated  to  the  point  where  a  sensor  is  located  should  be  taken  into  ac¬ 
count.  In  the  ordinary  art  gallery  problem,  these  effects  are  reduced  to  the 
geometric  shadow  zone:  any  point  within  the  geometrical  shadow  of  an  oc¬ 
clusion  for  a  certain  camera  is  invisible  for  this  camera.  In  probabilistic 
terms,  this  means  that  the  probability  of  detection  is  zero  for  the  shadow 
zone  and  is  one  for  the  field  of  view  before  the  occlusion.  This  idealized 
picture  is  a  consequence  of  the  geometrical  ray  approximation  for  signal 
propagation.  The  probabilistic  framework  is  capable  of  describing  more 
complicated  situations.  The  probability  of  detection  depends  on  the  signal 
and  noise  probability  density  functions  (Burdic  1984).  The  signal  and 
noise,  while  traveling  toward  a  sensor,  are  subjected  to  all  wave  phenom¬ 
ena,  such  as  scattering,  attenuation  (generally  inhomogeneous),  reflection 
from  objects  and  boundaries,  refraction,  interference  with  other  signals, 
and  diffraction.  Obstacles  and  terrain  variations  within  a  specified  area  of 
coverage  may  have  different  sizes,  shapes,  and  attenuation  factors  (corre¬ 
sponding  to  totally  obscured,  semitransparent,  and  transparent  areas).  As 
a  result,  the  probability  of  detection  may  be  a  highly  complicated,  anisot¬ 
ropic,  inhomogeneous  function  of  sensor  and  source  locations. 

The  aim  of  this  technical  report  is  twofold.  First,  two  multimodal  sensor 
placement  strategies  that  satisfy  the  specified  coverage  preferences  and 
take  into  account  the  terrain-based  and  signal-propagation  effects  are  pre¬ 
sented.  The  coverage  preferences  are  formulated  by  the  direct  assignment 
of  the  desired  probability  of  detection  to  each  point.  For  concreteness,  the 
two-dimensional  problem  is  considered  in  this  report,  although  the  devel¬ 
oped  approaches  can  be  applied  in  the  three-dimensional  space.  The  first 
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strategy,  called  Strategy  l,  is  aimed  at  finding  the  minimal  number  of  sen¬ 
sors;  it  is  heuristic  in  nature  and  places  sensors  one-by-one  in  the  worst 
covered  spots  (see  section  3  for  a  strict  definition  of  such  spots).  This 
strategy  is  similar  to  that  described  in  Dhillon  and  Chakrabarty  (2003) 
with  three  distinctions:  1)  the  probability  of  detection  in  the  current  study 
is  modeled  more  realistically,  taking  into  account  the  effects  of  signal 
propagation,  rather  than  assuming  the  exponentially  decaying  law,  2)  the 
coverage  preferences  are  incorporated  into  the  placement  strategy  rather 
than  used  as  a  stop  criterion  for  a  placement  algorithm,  and  3)  two  criteria 
are  used  simultaneously  to  decide  what  is  the  best  place  for  the  next 
placement.  Strategy  1  yields  a  rapid  but  suboptimal  solution  (that  is,  it  may 
not  find  the  minimal  number  of  sensors  required).  The  second  strategy, 
called  Strategy  2,  is  based  on  a  solution  of  the  strictly  posed  binary  linear 
programming  optimization  problem  (Sierksma  2002).  Such  a  formulation 
for  the  camera  placement  (a  nonprobabilistic  sensor  framework)  is  pre¬ 
sented  in  Erdem  and  Sclaroff  (2004).  This  approach  allows  one  to  mini¬ 
mize  total  costs  of  the  sensor  network  by  specifying  a  cost  function.  The 
cost  function  reflects  costs  associated  with  the  placement  of  a  particular 
sensor  to  a  particular  place.  It  may  be  the  actual  cost  of  a  sensor,  the  cost 
to  mount  and  operate  a  sensor  at  a  particular  location,  or  vulnerability  of  a 
sensor.  The  cost  function  also  could  be  chosen  to  minimize  the  total  num¬ 
ber  of  sensors  (that  is,  the  problem  will  be  reduced  to  the  same  problem 
Strategy  1  aims  to  solve).  Whatever  cost  function  is  chosen,  Strategy  2 
yields  the  global  optimal  solution,  that  is,  there  is  no  other  configuration  of 
sensors  that  would  satisfy  the  coverage  requirements  with  lesser  costs. 
However,  at  high  spatial  resolutions,  the  optimization  problem  becomes 
very  demanding  of  computational  time  (the  problem  belongs  to  the  nonde- 
terministic  polynomial,  more  precisely,  NP-complete,  class  of  complexity 
that  requires  significant  computational  resources;  Sierksma  [2002]).  In 
practice,  a  strict  solution  of  this  problem  at  large  matrix  dimensions  is  not 
obtainable.  For  example,  in  Erdem  and  Sclaroff  (2004),  the  set  of  possible 
camera  locations  was  quite  sparse  to  allow  one  to  solve  the  optimization 
problem.  In  this  technical  report,  a  fast  algorithm  for  finding  an  approxi¬ 
mate  solution  to  this  problem  in  the  case  of  large  dimensions  is  developed. 

Second,  a  safe-mode  concept  for  sensing  is  presented  in  this  study.  This 
concept  takes  into  account  that  the  probability  of  detection  of  a  given  sen¬ 
sor  depends  on  the  sensor's  operational  conditions,  such  as  weather,  ter¬ 
rain,  and  time  of  day.  The  same  operational  conditions  may  worsen  the 
probability  of  detection  of  one  sensor,  but  improve  others.  Defining  safe 
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modes  as  having  probability  of  detection  of  sensors  within  certain  ranges, 
one  can  calculate  the  probability  that  at  least  one  of  the  safe  modes  occurs 
under  various  operational  conditions.  A  more  detailed  explanation  of  this 
concept  is  given  in  section  6.  The  concept  is  an  adaptation  of  the  voting 
logic  inference,  presented  in  Klein  (2007),  to  the  probabilistic  sensor  per¬ 
formance  framework. 

The  report  is  organized  as  follows.  Section  2  describes  how  the  probability 
of  detection  can  be  modeled  to  take  into  account  signal  propagation  ef¬ 
fects.  Strategy  1  is  presented  in  section  3.  Section  4  is  devoted  to  Strategy  2 
and  the  algorithm  for  an  approximate  solution  of  the  binary  linear  pro¬ 
gramming  problem.  Section  5  presents  the  case  study  of  the  placement 
strategies  for  various  coverage  preferences  and  obstacles.  The  safe-mode 
concept  is  presented  in  section  6.  Conclusions  are  drawn  in  section  7. 
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2  Probability  of  Detection  with  Signal 
Propagation  Effects 

To  take  into  account  signal  propagation  effects,  the  following  physical 
model  is  adopted  in  this  technical  report.  An  object,  subjected  to  detection 
and  located  at  r0 ,  generates  multimodal  signatures,  called  signals  and  de¬ 
noted  s(rC),t) ,  t  being  the  time.  These  signals  have  random  amplitude, 
shape,  and  duration.  While  traveling  to  a  sensor  location,  the  signals  are 
subjected  to  scattering,  attenuation,  and  other  wave  propagation  effects; 
they  may  be  corrupted  by  noise  n(rR,t) ,  so  that  an  imperfect  signal 
u(rR,rQ,t)  =  s(rR , r 0,t)  +  n(rR , t )  is  received  by  a  sensor  located  at  .  The 

strict  solution  of  what  signal  will  be  received  at  a  sensor  location  would  in¬ 
volve  a  solution  of  the  wave  equation  with  properly  specified  boundary  and 
initial  conditions,  including  conditions  on  obstacles  and  terrain  inho¬ 
mogeneities,  which  is  beyond  the  scope  of  this  report.  Instead,  some  basic 
effects  of  signal  propagation  in  the  media  with  losses  are  taken  into  ac¬ 
count,  which,  nevertheless,  allow  one  to  capture  a  realistic  dependence  of 
the  probability  of  detection  on  distance.  Note  that  this  model  is  not  aimed 
at  enhancing  or  competing  with  a  strict  theory  of  signal  propagation  in 
random  media  (Tatarskii  1971;  Ishimaru  1978;  Ostashev  1997).  In  this  the¬ 
ory,  the  goal  is  to  find  relationships  between  statistical  properties  of  ran¬ 
dom  media  and  statistical  properties  of  the  propagated  signal.  In  the 
model  presented  in  this  technical  report,  the  medium  is  considered  to  be 
nonrandom  but  inhomogeneous;  that  is,  it  may  contain  obstacles  or  other 
objects  that  are  not  random  in  nature,  unlike  turbulence.  Taking  turbu¬ 
lence  into  account  could  be  a  further  enhancement  of  the  proposed  model. 

The  probability  of  detection  Pd(rR,r0)  of  a  signal  emitted  at  r0  by  a  sensor 
located  at  depends  on  the  energy  of  the  signal  and  noise  and  can  be  cal¬ 
culated  as  follows  (Burdic  1984): 


PArR’ro)  =  \  Pu(rR,ro)MdW’ 

r 


(i) 
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where  pv  is  the  probability  density  function  of  energy  U(rR,rQ)  of  the 
noisy  signal  u(rR,rQ,t),  and  y  is  the  detection  threshold  producing  the  pre¬ 
scribed  probability  of  false  alarm  Pfa  (rA, ) : 


Pfa(rR)  =  \PN(rR)(w)dw,  (2) 

r 


where  pN  is  the  probability  density  function  of  energy  N(rR )  of  the  noise 
n(rR,t) .  Although  Equation  2  allows  one  to  specify  Pfa  as  a  function  of  rR, 

a  usual  choice  would  be  to  have  a  constant  probability  of  false  alarm  for 
each  sensor  located  at  any  rR  within  an  area  being  protected.  To  find  y, 

one  may  notice  that  it  satisfies  a  requirement  that  the  cumulative  probabil- 

7 

ity  density  function  PM(ri;)(y)  =  j pNirR](w)dw  would  be  equal  to  l-P/a(rfi). 

—00 

Then: 


y~ ^/a(ri?))’  (3) 

where  FN\rR)  stands  for  the  inverse  cumulative  probability  density  func¬ 
tion. 

The  signal  s(rR,rQ,t)  can  be  presented  as  s(r^,ro,0  =  s(rR,rQ)s(t)  because 
the  wave  equation  allows  separation  of  variables  (in  the  absence  of  phase 
speed  scatterers  or,  in  practice,  when  they  are  negligibly  weak;  such  an  ap¬ 
proximation  is  used  in  x-ray  tomography,  for  example,  where  only  the  sig¬ 
nal's  attenuation  is  taken  into  account).  For  example,  for  a  plane  wave 
propagating  in  the  two-dimensional  isotropic  homogeneous  medium  with 
losses  (see,  e.g.,  Ostashev  [1997]): 

s(rfi ,  rG)  =  exp(-  a(oj)r  +  ik(a>)r)/ 4r  (4) 

and 

s(t)  =  Aexp(-zV^),  (5) 

where  r  is  a  distance  between  the  transmission  and  reception  r  =  |rfi  -rG|, 
a(co)  is  the  amplitude  attenuation  factor,  which,  in  general,  may  depend 
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on  the  frequency  co,  i  is  the  imaginary  unit,  k(co)  is  the  wave  number,  and 

A  is  the  amplitude.  This  example  assumes  a  harmonic  signal  with  a  single 
frequency  co.  For  more  complicated  signals,  s(t)  is  expressed  in  terms  of 

the  Fourier  series  or  integral.  Note,  however,  that  for  each  of  the  harmon¬ 
ics  in  the  Fourier  spectrum,  the  spatial  term  s(r^,r0)  can  be  presented  in 
the  form  given  by  Equation  4  with  possibly  distinct  a(co)  and  k(co).  As  a 
result,  the  signal's  energy  detected  by  a  sensor,  5"(rs,r0),  which  is  propor- 

tionalto  |s(rfi,r0,r)|  =  s(rfi,ro4>’(r*>ro>0,  where  *  denotes  complex  conjec¬ 
ture,  takes  the  form: 


S(rR  ,rQ)  =  At  |s(rfi ,  r0  )|2  J  |^(r)|2  dt  = 

0 

cxp(—  fir)  c 

—  ^0? 
r 


where 


T 

S0=4\\s(t)\2dt, 

0 


(6) 


(7) 


fi(co)  =  2  a(co)  is  the  energy  attenuation  factor,  4,2  is  a  coefficient  of  pro¬ 
portionality,  which  may  reflect  the  sensitivity  of  a  sensor  and  the  aperture 
of  the  signal  reception,  T  is  the  duration  of  the  signal,  and  S0  is  the  en¬ 
ergy  contained  in  the  temporal  part  of  the  signal  s(rR,rQ,t) .  Because  the 
duration  and  the  amplitude  of  the  signal  s(t)  are  random,  S0  is  a  random 

non-negative  variable.  Its  probability  distribution  function  is  not  easily  ob¬ 
tainable.  At  a  fixed  T  and  normally  distributed  s(t)  oc  .¥(/// /),cr((/)},  its 

probability  distribution  function— after  the  integral  is  digitized  and  ap¬ 
proximated  by  a  sum  with  m  terms— would  be  proportional  to  the  non¬ 
central  x1  distribution  with  m  degrees  of  freedom  and  the  parameter  of 

noncentrality  A  =  X™i  (/A,  / ‘A/  )2 ■  However,  the  authors  are  not  aware  of  a 
defined  probability  distribution  function  for  S0  with  random,  arbitrary 
distributed  m,  and  arbitrary  distributed  s(t).  In  this  technical  report,  it  is 
assumed  that  S0  can  be  well  described  by  the  normal  distribution 
N(/jSii ,  (jSii  )  with  positive  /us  >  3<t5  .  It  follows  from  Equation  6  that 
A ( ,  r0 )  is  also  normally  distributed  with  the  mean  jus  and  standard  de¬ 
viation  as  given  by  the  following  formulas: 
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Ms(rR>ro) 


exp(-  fk) 


vs(r. 


R  ’  r0  )  — 


exp(-  fir) 


<T, 


(8) 


Note  that  for  some  sensors  (e.g.,  seismic  or  magnetic)  the  power  of  the  r 
in  the  denominator  can  be  greater  than  one. 

If  energy  attenuation  is  spatially  nonuniform,  that  is  /?  =  J3(a>,  r),  one  can 
calculate  integral  attenuation  B(rR,rQ)  using  a  ray  approximation  and  tak¬ 
ing  the  integral  along  a  ray's  path: 

B(co,rR,rQ)=  J  j3(co,r)dl,  (9) 

LrorR 


where  L  is  the  length  of  a  ray  path  from  r0  to  rR.  Then,  this  B(co,rR,rQ) 
should  be  used  instead  of  fir  in  Equation  6.  Note  that  for  straight  paths 
and  uniform  ft,  Equation  9  yields  the  fk  factor. 

The  noise  energy  N(rR)  is  modeled  taking  into  account  the  following  con¬ 
siderations.  The  signal  is  corrupted  by  noise  n(rR,t)  at  the  time  and  point 
of  receiving.  The  energy  of  this  noise  accumulated  in  the  signal  of  the  du¬ 
ration  T  is: 


(10) 


where  the  constant  A\  is  a  coefficient  of  proportionality  and  may  reflect 

the  sensitivity  of  a  sensor  to  noise.  Following  the  same  arguments  as  for 
the  signal's  energy,  N( rR)  is  assumed  normally  distributed  with  the 
mathematical  expectation  juN(rR)  and  standard  deviation  a  N (rR),  such 
that  juN  >  3a N. 

Finally,  the  energy  of  the  noisy  signal  u(rR , r 0,t)  =  , rG,  t)  +  n(rR ,t)  is: 

T 

u(rR  T0 )  =  J  \u(rR ,  rQ ,  0|2  dt  =  ,  rQ )  +  N(rR ).  (11) 

0 


In  Equation  11,  the  integrals  over  cross-terms  s*n  and  sn  yield  zero  be¬ 
cause  the  signal  and  the  noise  are  assumed  to  be  entirely  incoherent.  This 
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fact  is  well  known,  for  example,  in  optics  where  the  intensities  of  two  inco¬ 
herent  sources  are  summed  algebraically  in  contrast  to  coherent  sources 
(Balanis  1989).  Because  S  and  N  are  independent  normally  distributed 
quantities,  U  is  also  normally  distributed  with  the  mean  pu  and  variance 

a\j  given  by  the  following  formulae: 


Mu (r« , r0 )  =  Ms (r« A0  )  +  Mn (r«  h  c4(r*,r0)  =  (J2s(rR,r0)  +  a2N(rR).  (12) 

At  large  distances  from  a  source,  when  r  — >•  00,  the  signal's  mean  and  vari¬ 
ance  equal  zero  (see  Equation  8),  so  that  pure  noise  is  detected  by  a  sen¬ 
sor.  In  this  case,  pL,  =  juN  and  av  =uN,  as  follows  from  Equation  12,  the 

probability  density  functions  pa  and  pN  coincide,  and  Pd  =  Pfa  (all  detec¬ 
tions  are  false  alarms).  Figure  1  shows  positions  of  three  probability  den¬ 
sity  functions,  for  the  signal,  noise,  and  noisy  signal,  when  a  sensor  is  not 
infinitely  distant  from  a  source. 


Figure  1.  Probability  density  functions  of  the  signal  energy  ps  (left),  noise  energy  pw  (middle), 
and  noisy  signal  energy  pu  (right).  The  horizontal  axis  represents  dimensionless  energy  w.  The 
source  is  quite  distant  from  a  receiver  so  that  its  energy  is  less  than  the  noise  energy.  The  y 
is  found  as  a  quantile  of  pn  corresponding  to  a  prescribed  probability  of  false  alarm.  The 
probability  of  detection  is  found  by  the  integration  of  pu  from  y  to  oo. 
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For  finite  r ,  the  amplitude  and  variance  of  a  signal  are  not  zeros  (see 
Equation  8),  juu>  juN  and  au>  aN,  which  causes  the  pLJ  curve  to  move  to 

the  right  relative  to  the  noise  distribution  function,  and  Pd  becomes 
greater  than  Pfa .  In  Figure  l,  w  (the  horizontal  axis)  denotes  dimen¬ 
sionless  energy  and  the  vertical  line  shows  the  value  of  y .  At  some  short 
distances  to  a  source,  the  pv  curve  lies  entirely  in  the  right-hand  side  area 
from  pN  and  Pd  =  1 . 


A  one-dimensional  spatial  section  of  the  probability  of  detection  of  a  signal 
emitted  at  x,  i.e.,  r0  =  (x,0) ,  by  a  sensor  located  at  rR  =  (0,0)  is  shown  in 

Figure  2. 


-2  -1.5  -1  -0.5  0  0.5  1  1.5  2 

X 


Figure  2.  One-dimensional  spatial  section  of  the  probability  of  detection  of  a  signal  emitted  at 
r0  =  (x,0)  by  a  sensor  located  in  the  origin. 

For  each  x  c=  [-2,2],  the  mean  and  variance  of  the  signal's  energy  were  cal¬ 
culated  in  accordance  with  Equation  8  with  r  =  |rx  -rG|  =  |x|  and  constant 
P .  These  mean  and  variance  define  the  location  and  shape  of  the  red  curve 
(the  signal)  in  Figure  l,  while  the  blue  curve  (the  noise)  does  not  move  (the 
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noise's  mean  and  variance  are  independent  of  the  distance  between  a 
source  and  a  receiver  by  assumption).  The  resulting  black  curve  (the  noisy 
signal)  also  moved  and  changed  its  shape  because  its  mean  and  variance 
changed  in  accordance  with  Equation  12.  The  y  was  determined  from  the 

blue  curve  (Equation  3)  at  probability  of  false  alarm  Pfa  =  10  6.  The  prob¬ 
ability  of  detection  (the  vertical  axis)  was  obtained  by  the  integration  of 
the  black  curve  from  this  y  to  infinity  (Equation  1).  As  one  can  see,  the 

spatial  dependence  of  the  probability  of  detection  falls  quite  suddenly  be¬ 
yond  a  certain  area  of  visibility  (the  detection  with  a  high  degree  of  prob¬ 
ability).  This  dependence  is  not  likely  due  to  the  exponential  or  Gaussian 
function.  A  numerical  experiment  with  the  noncentral  probability  den¬ 
sity  function  for  U( rR,rQ)  and  N(rR )  (at  a  fixed  signal  duration  T  )has 

been  carried  out  as  well.  The  probability  of  detection  obtained  was  similar 
to  that  shown  in  Figure  2  with  minor  quantitative  distinctions  while  the 
numerical  instability2  and  computational  time  were  significantly  greater. 
Thus,  the  Gaussian  probability  distributions  do  not  change  principal  re¬ 
sults  and,  therefore,  will  be  used  in  this  technical  report. 


2  Calculations  of  the  inverse  noncentral  chi-squared  distribution  functions  involve  a  summation  of  infi¬ 
nite  series,  which  have  poor  convergence  and  precision  at  certain  combinations  of  their  arguments 
(depending  on  the  chosen  probability  of  false  alarm  and  signal  and  noise  magnitudes  and  variances). 
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3  Strategy  1:  “Cover  the  Worst  First” 

The  strategy  described  in  this  section  is  similar  to  that  reported  in  Dhillon 
and  Chakrabarty  (2003).  It  belongs  to  the  class  of  "greedy"  algorithms, 
meaning  that  only  one  sensor  at  a  time  is  being  optimized.  Although  such 
algorithms  do  not  guarantee  that  the  final  sensor's  network  will  be  opti¬ 
mal,  they  are  widely  used  for  their  speed.  Strategy  1  is  heuristic  in  nature 
and  based  on  the  idea  that  sensors  are  placed  one-by-one,  each  in  the 
place  where  it  is  most  needed.  In  Dhillon  and  Chakrabarty  (2003),  this 
place  is  determined  as  the  place  with  the  minimal  probability  of  detection 
provided  by  a  current  sensor  network.  The  probability  of  detection  is  as¬ 
sumed  to  be  described  by  the  exponential  function  exp(-ar).  The  coverage 

preferences  are  formulated  as  the  desired  probability  of  detection  at  each 
spatial  point.  This  allows  one  to  assign  a  higher  desired  probability  of  de¬ 
tection  to  areas  of  high  importance  and  a  lower  one  to  areas  of  low  impor¬ 
tance.  If  there  are  no  specific  preferences,  one  can  assign  a  constant  de¬ 
sired  probability  of  detection  to  all  points.  If  there  is  a  certain  area  with 
higher  coverage  preferences,  then  one  continues  to  add  sensors  until  the 
preferences  are  satisfied.  In  other  words,  the  preferences  serve  as  a  stop 
criterion  and  are  not  incorporated  in  the  placement  strategy. 

Strategy  1  described  in  this  section  has  three  enhancements.  First,  the 
preferences  are  incorporated  into  the  strategy.  This  is  accomplished  by  as¬ 
suming  that  the  worst  covered  spot  is  that  where  the  discrepancy  between 
the  current  and  desired  coverage  is  the  largest.  Second,  a  more  realistic 
function  of  spatial  dependence  of  the  sensor's  probability  of  detection  is 
used.  This  function  (see  Figure  2)  falls  quite  suddenly  beyond  a  certain 
area  and,  therefore,  there  are  many  spatial  points  that  are  covered  equally 
poorly.  Therefore,  it  is  worthwhile  to  have  another  criterion  to  decide 
which  point  among  those  is  the  best  in  which  to  place  a  sensor.  A  geomet¬ 
rical  criterion  is  used  in  our  approach,  which  is  the  third  enhancement.  It 
is  assumed  that,  among  those  points,  the  most  distant  point  from  all  sen¬ 
sors  is  the  best.  If  there  are  many  best  points,  one  can  take  any,  for  exam¬ 
ple,  the  first  one.  Taking  these  considerations  into  account,  Strategy  1  can 
be  formulated  as  follows. 
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Step  l.  Calculate  the  current  probability  that  at  least  one  sensor  in  a  net¬ 
work  detects  a  source,  jP(r0;rA,,...,rA,  ),  where  M  is  the  current  number  of 

sensors,  (r*  ,...,rR  )  are  their  locations,  and  rQ  =(x,y)  indicates  the  possi¬ 
ble  locations  of  a  source  (that  is,  each  spatial  point  within  a  protected 
area).  P  is  considered  as  a  function  of  r0  while  rR  ,...,rR  are  fixed  pa¬ 
rameters.  To  simplify  notations,  these  parameters  will  be  omitted  when¬ 
ever  this  will  not  confuse  the  derivations.  It  is  assumed  that  sensors  detect 
a  source  independently.  This  assumption  allows  one  to  calculate 
Pc  (r0 ;  rR  , . . . ,  rR  )  very  efficiently  through  the  probability  of  misdetection 

Pmd  =  1  -  P  (not  to  be  confused  with  the  probability  of  false  alarm  Pfa ). 

Pmd(r0-rRi,...,rRu)  is  the  probability  that  all  M  sensors  fail  to  detect  a 

source  simultaneously.  Because  sensors  are  assumed  to  be  independent, 

Pmd  can  be  expressed  in  terms  of  the  probability  of  misdetection  of  a  single 

sensor  Pmd( r0;r*(),  i  = 

M 

Pmd  (r0;rRi,..„  vRm  )  =  Yl  Pmdi  (rc ;  ),  (13) 


where  the  probability  of  misdetection  of  a  single  sensor  is  calculated  using 
the  single  sensor's  probability  of  detection  Pd( rQ;rR  )  given  by  Equation  l: 

pmdi  (ro ;  ^ )  = 1  -  pd  (ro ;  rR,  )•  ( : 14 ) 

Then,  the  probability  that  at  least  one  sensor  of  M  would  detect  a  source 
can  be  found  easily: 


Pc  =  \-Pmd-  (15) 

If  no  sensors  are  placed  yet  ( M  -  0 ),  one  sets  P  =  0  for  each  point. 

Step  2.  Calculate  the  discrepancy  s  between  the  current  and  desired 
probabilities  of  detection: 

£(r0)  =  Pc(ro)-Ppr(r0),  (16) 

where  Ppr{ r0)  is  the  desired  probability  of  detection  of  a  source  at  rQ.  In 

principal,  this  criterion  alone  already  allows  one  to  find  the  best  location 
for  another  sensor  ,  which  would  correspond  to  the  minimum  of  s : 
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where  Ppr( rQ)  is  the  desired  probability  of  detection  of  a  source  at  r0.  In 

principal,  this  criterion  alone  already  allows  one  to  find  the  best  location 
for  another  sensor  rR  ,  which  would  correspond  to  the  minimum  of  a : 

rflu+l  =  argmin[£-(r0)].  (17) 

If  there  are  many  points  where  a  reaches  its  minimum,  one  takes  the 
point  in  the  middle  of  that  set.  This  approach  is  used  in  the  numerical  ex¬ 
periment,  described  in  this  report,  in  the  absence  of  other  sensors  (that  is, 
to  find  the  best  location  for  the  first  sensor).  If  there  are  other  sensors 
(M  >  0 ),  criterion  a  can  be  enhanced  by  a  geometrical  criterion,  as  shown 
in  Step  3. 

Step  3.  Calculate  criterion  d0.  Criterion  <70  takes  into  account  current  lo¬ 
cations  of  sensors  and  defines  the  best  point  for  another  sensor  as  the 
most  remote  point  from  all  existing  sensors.  More  precisely,  for  each  spa¬ 
tial  point  r0  of  possible  source  locations,  one  finds  the  distance  to  the 
nearest  sensor  dQ(r0)  =  min {d(rQ,rR  ),..., d(r0,rR  )},  where  d(rQ,rR)  are  dis¬ 
tances  from  rQ  to  sensors  located  at  .  Then,  according  to  criterion  d0, 
the  best  point  is  that  where  d0  reaches  its  maximum: 

rRM+l  =  argmax[d(r0)].  (18) 

Criterion  d0  alone,  as  well  as  criterion  a,  can  be  used  to  place  sensors.  In 

this  case,  the  locations  can  be  found  in  advance  purely  from  geometrical 
considerations,  without  any  characteristics  of  sensors  required. 

Step  4.  Calculate  a  supercriterion  C  that  combines  both  the  a  and  d0  cri¬ 
teria  together.  To  do  this,  one  needs  to  normalize  a  and  d0  so  that  they 

have  the  same  range  of  variability.  The  following  normalization  yields  the 
range  [0,1]  for  normalized  a  and  d0  (C,  and  C2,  correspondingly): 

max[g(r0)]-g(r0) 
max[4r0)]-min[4r0)]’ 

C2(rQ)=  do(r°}  , 
max[d0(ro)] 


(19) 


ERDC/CRREL  TR-08-24 


15 


Supercriterion  C  is  nothing  but  the  sum  of  C,  and  Cv  Thus: 

r*„+1  =  ar8  max[C2  (r0  )  +  C,  (rG )].  (20) 

Supercriterion  C  operates  as  follows.  If,  for  example,  there  are  two  points 
with  the  same  discrepancy  between  the  current  and  desired  probability  of 
detection,  then  the  most  distant  point  from  all  sensors  will  be  selected. 

And  vice  versa,  if  there  are  two  points  equally  remote  from  all  other  sen¬ 
sors,  the  one  with  the  minimal  e  (or,  equivalently,  with  the  maximal  C,, 

which  corresponds  to  the  higher  discrepancy  between  the  current  and  de¬ 
sired  probabilities)  will  be  selected.  If  there  are  many  points  with  the  same 
value  of  C,  any  could  be  selected,  for  example,  the  first  one.  To  vary  the 
importance  of  C2  or  Cl  in  the  final  decision,  one  can  assign  different 
weight  factors  W2  and  Wl  to  criteria  C2  and  C, . 

Step  5.  Once  the  new  placement  point  is  known,  the  last  step  is  to 

determine  what  type  of  sensors  should  be  placed  there.  Because  the  goal  is 
to  find  the  minimal  number  of  sensors,  then  a  sensor  with  a  wider  area  of 
coverage  (if  placed  at  )  will  be  the  choice.  Now  the  current  number  of 

sensors  is  M  + 1,  and  one  repeats  the  algorithm  from  Step  1  until  all  points 
are  covered  with  the  desired  probability  of  detection  (or  greater)  or  the 
number  of  available  sensors  is  reached. 

The  described  strategy  has  several  advantages.  First,  it  is  very  fast.  On  a 
spatial  grid  81x81  points,  it  takes  a  fraction  of  a  second  to  place  a  virtually 
unlimited  number  of  sensors  (on  Intel  Xeon  2.66  GHz  PC  with  4  Gb  of 
RAM).  Second,  to  place  a  new  sensor,  one  does  not  need  to  remount  and 
replace  already  existing  sensors.  In  other  words,  the  strategy  allows  an  ex¬ 
pansion  of  an  existing  network.  If,  for  example,  a  current  area  being  cov¬ 
ered  becomes  larger,  one  can  take  an  existing  sensor  network  as  a  starting 
point  and  expand  it  to  match  the  new  area.  And  third,  the  strategy  can  be 
limited  by  a  maximal  number  of  available  sensors  even  if  the  desired  prob¬ 
ability  of  detection  was  not  reached.  As  one  will  see  below,  the  strict  solu¬ 
tion  does  not  allow  such  a  limitation. 

The  disadvantage  of  this  strategy  is  that  it  is  suboptimal.  Although  the  so¬ 
lution  found  is  much  better  than  a  random  placement,  as  demonstrated  in 
Dhillon  and  Chakrabarty  (2003)  for  the  two-dimensional  case,  it  may  not 
have  the  widest  possible  area  of  coverage  for  a  given  number  of  sensors; 
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or,  if  there  is  no  limitation  on  the  total  number  of  sensors,  the  solution 
may  not  yield  the  minimal  possible  number  of  sensors  to  satisfy  the  cover¬ 
age  preferences.  To  illustrate  this  point,  imagine  that  one  needs  to  cover 
uniformly  a  one-dimensional  segment  by  a  minimal  number  of  sensors. 
Suppose  the  sensors  are  identical  and  have  the  radius  of  visibility  equal  to 
one-third  of  the  segment's  length.  According  to  Strategy  l,  the  first  sensor 
will  be  placed  in  the  center  of  this  segment,  covering  two-thirds  of  its 
length,  the  second  sensor  will  be  placed  in  one  of  the  segment's  vertices 
(because  they  are  most  remote  from  the  center  and  have  the  worst  cover¬ 
age),  and  the  third  one  will  be  placed  in  the  other  vertex.  Thus,  the  answer 
would  be  that  three  sensors  are  required.  Obviously,  the  strict  optimiza¬ 
tion  would  place  only  two  sensors  (for  example,  in  the  positions  that  divide 
the  segment  into  three  equal  parts). 
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4  Strategy  2:  “Do  the  Job  with  the 
Minimal  Cost” 

Strategy  2  was  developed  to  overcome  the  nonoptimality  of  Strategy  1.  It  is 
based  on  a  solution  of  the  strictly  posed  optimization  problem,  and  guar¬ 
antees  determination  of  a  global  minimum  (if  it  exists)  of  costs  specified 
by  a  cost  function.  The  cost  function  may  reflect  the  actual  costs  of  sen¬ 
sors,  or  other  disincentives,  such  as  the  total  number  of  sensors,  the  sen¬ 
sors'  vulnerability,  or  their  reliability.  For  fine  spatial  resolution,  which 
leads  to  large  matrix  dimensions,  however,  the  strict  solution  of  this  prob¬ 
lem  is  not  practically  obtainable.  The  same  strict  problem  formulation  for 
camera  placement  can  be  found,  for  example,  in  Erdem  and  Sclaroff 
(2004).  In  this  work,  the  authors  were  forced  to  have  a  rather  coarse  spa¬ 
tial  grid  for  possible  camera  locations  to  be  able  to  solve  the  problem.  In 
this  technical  report,  an  approximate  solution  of  the  strictly  posed  optimi¬ 
zation  problem  is  found  for  large  dimensions,  allowing  fine  resolution  for 
both  the  sensor  placement  and  the  area  to  be  covered. 

Strict  formulation  of  the  optimal  coverage  problem 

In  this  subsection,  it  is  shown  how  the  problem  of  optimal  coverage  in  the 
probabilistic  sensor  performance  framework  can  be  strictly  formulated  as 
the  binary  linear  programming  problem.  Let  us  introduce  a  column-vector 
p.  The  length  of  vector  p  equals  the  number  of  possible  sensor  locations 

K\  only  a  few  of  them  will  be  chosen  as  optimal  ones.  Furthermore,  let  us 

require  that  the  elements  of  this  vector  could  take  either  0  or  1  values 
only.  Then,  vector  p  can  be  treated  as  the  indicative  vector:  its  0  element 

in  the  k  th  position  shows  that  no  sensors  should  be  placed  in  the  k  th  pos¬ 
sible  point  of  location  while  1  in  this  position  indicates  that  a  sensor 
should  be  placed  there.  If  the  goal  is  to  find  a  minimal  number  of  sensors, 
then  the  optimal  vector  p0  would  have  a  minimal  number  of  ones.  Their 

positions  in  this  vector  indicate  where  to  place  sensors.  Mathematically, 
this  requirement  is  formulated  as  finding  a  minimum  of  a  scalar  linear 
function  F  =  f  rp,  where  f  is  a  column  vector  of  costs  associated  with  a 

particular  sensor  and  T  indicates  the  transposition.  For  a  minimal  num¬ 
ber  of  sensors,  f  =  [1;1;...;1],  where  the  semicolon  between  elements  de¬ 
notes  that  these  elements  are  arranged  in  a  column.  But  f  can  consist  of 
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other  values,  for  example,  fk  can  be  equal  to  the  actual  costs  of  placing  a 
sensor  in  the  k  th  possible  point  or  it  could  be  the  probability  that  a  sensor 
will  be  found  and  disabled  at  this  point.  From  this  perspective,  the  F 
function,  which  is  subject  to  minimization,  is  nothing  but  the  total  cost  of 
placing  the  sensors.  If  there  are  several  types  of  sensors,  then  there  is  a 
vector  of  possible  locations  for  each  of  them,  and  the  length  of  the  p  vector 

changes  accordingly.  For  example,  if  the  possible  locations  for  two  types  of 

sensors  are  the  same  (e.g.,  length  K ),  then  the  P  vector  would  be  length 
2 K,  and  the  number  and  positions  of  ones  would  indicate  how  many  and 
what  type  of  sensors  are  required.  The  cost  function  f ,  in  this  case,  also  is 
length  2 K  and  reflects  the  costs  for  the  first  and  second  types  of  sensors. 

To  impose  the  constraints  of  coverage,  let  us  note  that  requiring  the  prob¬ 
ability  of  detection  Pd  greater  or  equal  to  the  given  preferences  Ppr  is 

equivalent  to  requiring  the  probability  of  misdetection  Pmd  less  or  equal  to 
the  allowed  probability  of  misdetection:  Pmd  <  1  -  Ppr.  According  to  Equa¬ 
tion  13,  Pmd  is  a  product  of  the  individual  Pmd 's.  Therefore,  taking  the  loga¬ 
rithm  from  both  sides  of  this  inequality,  one  has  for  each  point  rQ : 


M 

X  ln  P’nd,  (rO  >  r«,  )  ^  ln0  -  Ppr  (rO  )}  (21) 

i= 1 

where  Pmd  (r0,rR  )  are  given  by  Equation  14.  In  this  context,  M  equals  the 
number  of  ones  in  the  vector  p.  If  the  total  number  of  possible  source  lo¬ 
cations  r0  is  Q,  the  last  equation  can  be  rewritten,  in  matrix  notations, 

Dp<b,  (22) 

where  b  =  ln(l-Ppr(r0)),  the  k  th  column  of  matrix  D  (size  Q*K)  repre¬ 
sents  the  logarithm  of  the  probability  of  misdetection  for  a  sensor  placed 
in  the  k th  point,  In Pmd  (rQ,rR ),  and  the  inequality  reads  as  an  element -by¬ 
element  comparison.  Note  that  the  number  of  terms  actually  selected  for 
the  summation  in  Equation  22  is  defined  by  the  number  and  position  of 
ones  in  the  vector  p. 

Thus,  the  problem  of  the  optimal  coverage  can  be  formulated  as  the  mini¬ 
mization  of  total  costs  with  constraints  that  reflect  the  coverage  prefer- 
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ences  and  restrict  p  to  be  a  binary  vector.  That  is,  one  seeks  an  optimal 
vector  p0  such  that: 


and 


p0  =argminfrp, 


Dp  <  b,  pk  =  {0,1},  k  = 


(23) 


(24) 


This  is  one  of  the  forms  of  the  binary  linear  programming  problem.  If  the 
problem  is  feasible  (that  is,  the  constraints  permit  at  least  one  solution) 
then  its  solution  is  guaranteed  to  be  the  global  minimum  of  the  total  costs 
(because  F  is  a  linear  function  of  p  ;  Sierksma  [2002]).  This  means  there 

is  no  other  configuration  of  sensors  satisfying  the  coverage  preferences 
that  would  decrease  the  total  costs  further.  In  the  example  with  a  one¬ 
dimensional  segment,  introduced  in  section  3,  a  strict  solution  of  this 
problem  would  place  only  two  sensors  within  the  segment. 

However,  for  large  dimensions  of  vector  p  and  matrix  D,  a  strict  solution 
is  difficult  to  obtain.  The  binary  linear  programming  problem  does  not 
have  a  quick  solution,  unlike  its  nonbinary  counterpart,  which  can  be 
solved  effectively  by  simplex  or  interior  point  methods.  To  give  insight  why 
the  binary  requirement  complicates  the  solution,  let  us  consider  a  basic 
idea  of  the  boundary  and  branch  method,  one  of  the  most  widely  used 
methods  for  the  binary  optimization  problems.  The  details  can  be  found, 
e.g.,  in  Sierksma  (2002).  At  first,  a  nonbinary  problem  is  solved  with  the 
constraints  that  all  elements  of  p  must  be  within  the  interval  [0, 1] .  Then, 

each  noninteger  element  of  this  solution  is  replaced  by  0 ,  which  generates 
one  branch  of  the  possible  binary  solution,  and  1 ,  which  generates  another 
branch  of  the  possible  binary  solution.  Therefore,  if  the  number  of  nonin¬ 
teger  elements  is  L,  the  total  number  of  branches  to  investigate  is  2L.  In 

this  method,  a  smart  way  to  investigate  and  cut  unnecessary  branches  is 
used,  based  on  the  constraints  (matrix  D) .  In  the  numerical  simulation 

described  in  this  report,  the  area  in  consideration  was  partitioned  into 
81x81  spatial  cells,  so  that  the  number  of  possible  sensor  locations  and 
possible  signal  emissions  was  812  =  6561.  Hence,  the  size  of  D  was  [6561  x 
6561]  and  the  length  of  p  was  6561.  In  the  worst  case,  all  elements  of  p, 

after  solving  the  nonbinary  problem,  can  be  noninteger,  and  one  will  need 
to  investigate  2  6561  branches.  Regardless  of  how  smart  the  algorithm  of 
branch  selection  is,  the  total  numbers  of  branches  and  constraints  are 
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large  enough  to  avoid  finding  any  solution  with  reasonable  computational 
effort.  An  algorithm  proposed  in  the  next  subsection  finds  an  approximate 
(almost  optimal)  solution  to  this  problem  in  less  than  one  minute  using  an 
Intel  Xeon  2.66  GHz  PC  with  4  Gb  of  RAM. 

Approximate  solution  for  minimal  number  of  sensors 

In  this  subsection,  an  algorithm  for  approximate  solution  to  the  binary  lin¬ 
ear  programming  problem  is  presented  for  the  case  when  f  =  [1;  1; . . . ;  1] ; 

that  is,  when  the  goal  is  to  minimize  the  total  number  of  sensors.  A  gener¬ 
alization  of  this  algorithm  for  other  f 's  will  be  given  in  the  next  subsec¬ 
tion.  The  approximate  solution,  presented  in  this  subsection,  converts 
Strategy  2  into  a  "greedy"  algorithm  (only  one  sensor  at  a  time  is  opti¬ 
mized).  However,  Strategy  2  places  sensors  in  a  significantly  different  way 
than  Strategy  1  as  will  be  demonstrated  in  Section  5. 

One  can  note  that,  for  the  binary  vector  p,  the  minimal  and  maximal  val¬ 
ues  of  total  costs  F  =  f  p  are  known  in  advance.  These  are  0  (when  p  =  0 ) 
and  K  (when  p  =  1 ).  The  only  reason  why  p  =  0  cannot  be  declared  as  the 
solution  is  the  constraint  Dp<b,  which  is  not  satisfied  by  this  p  (if  the 
constraint  is  satisfied,  then  the  optimal  solution  is  p0  =  0 ).  Note  that  if  this 
constraint  is  not  satisfied  by  p  =  1  either,  then  the  problem  is  infeasible, 
i.e.,  there  is  no  solution  to  this  problem  in  the  binary  space  for  p  .3  Thus, 
one  can  verily  very  quickly  whether  the  problem  is  feasible. 

If  the  problem  is  feasible,  one  needs  to  distribute  as  few  ones  as  possible 
among  positions  in  the  p  vector  (which  is  filled  in  by  zeros  otherwise). 

This  would  yield  a  strict  solution.  The  idea  of  the  algorithm  for  approxi¬ 
mate  solution  is  to  do  it  consecutively  to  avoid  investigating  all  possible 
combinations.  The  placement  of  a  single  1  in  any  position  in  p  is  equiva¬ 
lently  poor  from  the  standpoint  of  the  total  cost  F  because  it  will  equally 
increase  its  value  (from  0  to  1).  However,  the  constraints  will  be  satisfied 
differently.  Let  p  be  a  trial  solution  that  equals  p  with  the  extra  1  added. 

The  best  position  of  1  is  that  where  the  sum  over  the  positive  elements  of 
A  =  Dp  -  b  is  minimal: 


3This  means  the  coverage  preferences  are  not  satisfied  even  if  one  places  sensors  at  all  possible  points. 
Such  a  situation  is  likely  to  happen  when  K  «  Q,  i.e.,  when  a  large  area  is  supposed  to  be  covered 
with  a  few  allowed  sensor  locations. 
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(25) 

A+(p)=  (A9(p):  A(?(p)>0},  q  =  \,...,Q, 

where  /  is  the  number  of  points  where  the  coverage  is  not  satisfied.  If  at 
some  position  of  1  there  are  no  such  points  ( /  =  0 ),  then  all  coverage  pref¬ 
erences  are  satisfied,  and  the  optimal  vector  p0  equals  the  trial  vector  p 

where  this  happens.  Thus,  the  algorithm  for  an  approximate  solution  can 
be  formulated  as  follows. 

Step  l.  Verily  that  the  problem  is  feasible.  If  it  is  not,  then  there  is  no  so¬ 
lution.  If  it  is,  set  p0  =  0. 

Step  2.  Verify  that  there  is  at  least  one  point  where  the  coverage  is  not 
satisfied.  If  the  coverage  preferences  are  satisfied  everywhere,  return  p0.  If 

not,  proceed  as  follows.  Set  a  trial  vector  p  =  p0 .  Let  the  possible  number 

of  sensor  locations  be  K.  (These  possible  sensor  locations  correspond  to 
zeros  in  p).  Then,  consecutively  set  zero  elements  in  p  to  one  and  calcu¬ 
late  A  =  Dp  -  b.  Find  the  positive  elements  of  A  and  sum  them  over. 
Memorize  this  sum  in  E(k),  k  =  \,...,K.  If  at  some  k  there  will  not  be  posi¬ 
tive  elements  of  the  residual  ( /  =  0 ),  set  p0k  =  1  and  return  p0.  The  solu¬ 
tion  is  found.  If  not,  go  to  Step  3. 

Step  3.  Find  the  best  k0  where  E(k)  reaches  its  minimum: 

K  =argmin E(k).  Set  p0ko  =1. 

Step  4.  Exclude  unnecessary  points  from  possible  sensor  locations.  Obvi¬ 
ously,  at  least  point  k0  should  be  excluded  (one  does  not  need  to  place  two 

sensors  at  the  same  point).  But  it  is  not  the  only  possibility.  One  can  find 
all  spatial  points  where  the  coverage  requirements  are  already  satisfied  (by 
placement  of  one  sensor  at  a  certain  point  corresponding  to  k0).  Namely, 
these  points  are  defined  by  the  row  indices  of  matrix  D  at  which  Dp0  <  b. 

If  some  spatial  points  of  possible  sensor  locations  also  belong  to  the  area 
where  the  coverage  is  already  satisfactory,  then  all  these  points  should  be 
excluded  from  further  consideration  because  it  does  not  make  much  sense 
to  place  another  sensor  in  this  area.  This  significantly  reduces  the  number 
of  possible  sensor  locations  K  (note  that  the  spatial  point  corresponding 
to  k0  will  be  excluded  automatically).  Go  to  Step  2. 
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In  the  worst  case,  the  algorithm  will  converge  to  the  most  nonoptimal  so¬ 
lution  with  p0  =  1. 

This  strategy  can  be  viewed  as  walking  from  the  ultimate  optimal  point 
p  =  0  (where  the  cost  F  is  minimal  but  constraints  are  not  satisfied)  to  the 
most  nonoptimal  point  p  =  1  (where  the  cost  F  is  maximal  but  the  con¬ 
straints  are  guaranteed  to  be  satisfied)  along  the  ribs  of  a  K  -dimensional 
hypercube,  hoping  to  satisfy  the  coverage  requirements  before  the  most 
nonoptimal  point  is  reached,  as  illustrated  in  Figure  3. 


Figure  3.  Strategy  2  guides  a  walker  from  the  most  optimal  node  (0,0,0),  where  the 
constraints  are  not  satisfied,  to  the  most  nonoptimal  node  (1,1,1),  where  the  constraints  are 
satisfied,  without  jumping  to  nonadjacent  nodes.  The  colors  of  nodes  reflect  the  value  of  total 

cost  F. 

In  this  example,  K  =  3,  so  that  p  =  (pl,p2,p3).  The  color  of  the  nodes  cor¬ 
responds  to  the  value  of  the  total  cost  F.  The  lowest  cost  ( F  =  0 )  corre¬ 
sponds  to  the  blue  node  at  (0, 0, 0).  The  highest  cost  ( F  -  3 )  is  at  the  red 
node  at  (1,1,1).  One  needs  to  walk  from  (0,0,0)  to  (1,1,1)  without  jumping 
(this  is  the  requirement  of  the  consecutive  placement;  if  jumping  to  non¬ 
adjacent  nodes  were  allowed,  a  strict  solution  would  result).  Because  all 
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neighbors  around  the  blue  node  equally  increase  the  cost  F,  all  of  them 
are  of  the  same  color.  Among  these  equally  poor  directions,  the  algorithm 
selects  the  one  that  brings  a  walker  closer  to  satisfaction  of  the  problem 
constraints. 

This  simple  strategy  provides  surprisingly  good  results,  as  demonstrated 
in  section  5.  Among  many  possible  sensor  locations,  it  quickly  finds  a 
small  number  satisfying  the  coverage  preferences.  But  the  solution  may  be 
suboptimal,  as  a  consequence  of  the  consecutive  placement  (i.e.,  sensors 
are  placed  one-by-one).  In  the  one-dimensional  example,  introduced  in 
section  3,  the  first  sensor  would  be  placed  in  the  center  of  a  segment,  the 
second  one  in  the  uncovered  segment  of  one  of  the  halves,  and  the  third 
one  in  the  uncovered  segment  of  the  other  half.  At  first  glance,  this  is  very 
similar  to  the  results  of  Strategy  1,  but  the  distinction  will  be  clear  in  the 
two-dimensional  case,  as  demonstrated  in  section  5.  The  advantage  of  the 
consecutive  placement  is  that  the  total  number  of  sensors  available  can  be 
incorporated  into  the  stop  criterion  of  the  algorithm. 

Generalization  to  entirely  positive  or  negative  cost  functions 

Strategy  2  can  be  generalized  for  arbitrary  cost  functions  with  all  positive 
or  negative  elements.  This  is  a  very  practical  situation  because  the  costs  of 
different  sensors  (whatever  they  reflect)  are  usually  of  the  same  sign.  If  all 
elements  of  f  are  negative,  all  inference  of  Strategy  2  remains  the  same 
with  the  only  distinction  that  p  =  1  is  the  ultimate  optimal  solution,  p  =  0 

is  the  most  nonoptimal  one,  and  one  fills  the  ultimate  optimal  solution 
with  zeros  rather  than  with  ones.  In  terms  of  Figure  2,  a  walker  now  goes 
in  the  opposite  direction  from  (1,1,1)  to  (0,0,0).  For  concreteness,  the  posi¬ 
tive  cost  functions  will  be  considered  in  this  subsection:  f  >  0. 

In  this  case,  as  before,  p  =  0  and  p  =  1  are  the  ultimate  optimal  and 
nonoptimal  solutions,  correspondingly.  However,  different  trial  solutions 
in  Step  2  are  not  equivalent  now  from  the  total  cost  point  of  view  (all 
neighbor  nodes  would  have  different  colors  in  Figure  3).  Therefore,  one 
should  trade  off  between  growing  costs  and  coverage  satisfaction.  Step  2 
should  change  as  follows.  In  addition  to  E(k),  which  reflects  the  quality  of 

coverage,  one  needs  to  calculate  and  memorize  the  total  costs  F(k )  =  f  rp 
corresponding  to  the  trial  solution  p.  The  most  appealing  point  for  a  sen¬ 
sor  placement  would  be  that  where  the  coverage  is  maximal  while  the  cost 
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increase  is  minimal.  To  combine  two  criteria  together,  a  supercriterion  G 
can  be  introduced,  similarly  to  the  supercriterion  C  in  Strategy  l: 


G(k)  = 


E(k)-rmn[E(k)\ 
max[ii  (A:)]  -  mm\E{k)\ 


+ 


F (k)  -  mm[F (k)] 
max[F(k)\  -  minfEXA:)] 


(26) 


where  normalization  is  needed  to  make  these  two  competitive  factors 
equally  important.  To  make  them  unequally  important,  one  can  introduce 
weight  factors  for  each  of  the  terms  in  Equation  26.  The  best  position  is 
defined  now  as  k0  =  argminG(A:). 
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5  Numerical  Experiment 

The  goal  of  the  numerical  experiment  was  to  compare  Strategy  l  and 
Strategy  2  under  various  coverage  preferences  and  obstacle  locations.  The 
area  being  protected  (a  2-X-2  square  in  conditional  units)  was  partitioned 
into  81-X-81  grid  cells.  For  Strategy  1,  there  were  three  types  of  sensors, 
each  with  a  specific  attenuation  constant  /?  that  defines  the  radius  of  the 

sensor's  coverage  (see  Equation  8).  Moreover,  it  was  assumed  that  these 
/3 's  depend  on  a  placement  point.  For  example,  a  type  1  sensor  placed  at 

one  point  could  be  better  than  a  type  2  sensor  at  the  same  point,  although 
the  situation  may  reverse  for  another  point.  The  values  of  these  /? 's  were 
randomly  chosen  from  the  interval  [0,1].  A  sensor  of  type  1  had  a  denomi¬ 
nator  proportional  to  r 2  rather  than  r,  in  contrast  to  other  types.  The  en¬ 
ergies  of  signal  and  noise  were  modeled  in  accordance  with  section  2.  For 
Strategy  2,  only  one  type  of  sensors  was  available  with  random,  point- 
independent  j8.  This  allowed  one  to  better  illustrate  the  main  features  of 

Strategy  2  (the  features  of  Strategy  1  can  be  seen  without  this  restriction). 
As  mentioned  in  section  4,  Strategy  2  allows  a  strict  expansion  to  multi¬ 
modal  sensors.  The  goal  was  to  minimize  the  total  number  of  sensors, 
which  was  restricted  to  nine,  again  for  illustrative  purposes.  Without  this 
restriction,  both  strategies  would  eventually  satisfy  the  coverage  prefer¬ 
ences  (with  a  different  number  of  sensors).  The  probability  of  false  alarm 
Pfa  was  equal  to  10  6  for  each  sensor. 

Figure  4,  (a)  and  (b),  shows  the  locations  obtained  with  the  use  of  Strategy 
1  and  Strategy  2,  correspondingly,  for  the  case  of  empty  space.  The  color 
represents  the  probability  that  at  least  one  sensor  will  detect  a  source.  Dif¬ 
ferent  markers  in  Figure  4  (a)  correspond  to  different  types  of  sensors  se¬ 
lected  by  Strategy  1.  The  coverage  preferences  were  0.95  for  each  spatial 
point.  As  one  can  see,  neither  of  the  strategies  satisfied  the  coverage  pref¬ 
erences  by  the  means  of  nine  sensors.  However,  the  area  covered  by  Strat¬ 
egy  2  is  remarkably  larger  than  that  covered  by  Strategy  1.  This  can  be  in¬ 
terpreted  to  mean  that  Strategy  2  yields  more  optimal  results  than 
Strategy  1.  This  tendency  holds  in  more  complicated  cases,  as  seen  below. 

Figure  5  presents  the  sensor  locations  and  coverage  with  two  high-value 
areas  indicated  as  green  squares.  The  coverage  preferences  were  0.95  in¬ 
side  these  areas  and  0.8  outside.  Again,  Strategy  2  places  sensors  more 
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Figure  4.  Probability  of  detection  by  at  least  one  of  the  sensors.  Sensor  locations  obtained  by 
(a)  Strategy  1  and  (b)  Strategy  2.  The  coverage  preferences  were  0.95  for  each  point,  and  the 
total  number  of  available  sensors  was  limited  to  nine.  Different  markers  (a)  correspond  to 

different  sensor  types. 


Figure  5.  Coverage  provided  by  (a)  Strategy  1  and  (b)  Strategy  2.  The  green  squares  indicate 
high-value  objects.  The  coverage  preferences  are  0.95  inside  the  objects  and  0.8  outside. 

appropriately,  without  extensive  uncovered  spots,  as  one  can  see  in  Figure 

5(b). 

In  addition  to  high-value  areas,  there  may  exist  forbidden  areas  or  obsta¬ 
cles,  where  no  sensors  are  allowed.  One  can  distinguish  two  types  of  obsta¬ 
cles.  The  first  consists  of  obstacles  that  should  be  covered  with  certain 
preferences,  although  no  sensors  are  allowed  inside  obstacles.  The  second 
consists  of  obstacles  within  which  coverage  is  irrelevant.  One  can  think 
about  such  obstacles  as  if  they  have  the  coverage  preferences  equal  to  zero. 
In  this  case,  the  main  goal  of  coverage  will  be  satisfying  the  preferences 
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outside  the  obstacles.4  Figure  6  presents  the  sensor  placement  by  Strategy 
l  and  Strategy  2  in  the  presence  of  high-value  areas  (green  squares)  and 
obstacles  (black  squares).  Obstacles  of  the  first  type  can  be  seen  in  Figure 
6,  (a)  and  (b).  For  this  case,  the  coverage  preferences  were  0.95  inside  the 
green  squares  and  0.8  everywhere  else  including  the  interior  of  the  obsta¬ 
cles.  The  placement  for  the  obstacles  of  the  second  type  is  shown  in  Figure 
6,  (c)  and  (d).  For  this  case,  the  required  coverage  was  0.95  inside  the 
high-value  areas  and,  effectively,  zero  inside  the  obstacles. 


Figure  6.  Sensor  locations  in  the  presence  of  high-value  objects  (green  squares)  and 
obstacles  (black  squares),  (a)  Strategy  1,  obstacles  of  the  first  type,  (b)  Strategy  2,  obstacles 
of  the  first  type,  (c)  Strategy  1,  obstacles  of  the  second  type,  (d)  Strategy  2,  obstacles  of  the 

second  type. 


4ln  practice,  such  objects  may  represent  natural  nooks  or  objects  with  independent  (internal) 
guarding. 
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The  most  complicated,  yet  still  practically  important,  case  is  when  the  ob¬ 
stacles  coincide  with  the  areas  of  high  importance.  That  is,  one  would  like 
most  to  cover  a  certain  area  where  no  sensors  are  allowed.  Strategy  l  can¬ 
not  handle  this  case  because  this  strategy  places  a  sensor  into  the  worst 
covered  place,  which,  in  this  case,  happens  to  be  forbidden.  Therefore, 
Strategy  l  will  place  sensors  in  other  worst  covered  places,  irrelevant  to  the 
areas  with  high  importance.  The  only  way  to  make  Strategy  l  succeed  in 
this  case  is  to  specify  a  high-value  area  wider  than  the  forbidden  one,  so 
that  at  least  the  edge  would  be  accessible,  as  shown  in  Figure  7  (a).  In  con¬ 
trast,  Strategy  2  succeeds  without  this  additional  restriction.  It  can  handle 
the  situations  when  the  forbidden  areas  are  wider,  narrower,  or  exactly 
equal  to  the  high-value  areas.  For  example,  in  Figure  7  (b),  the  obstacle 
exactly  coincides  with  the  high-value  area  so  that  it  is  seen  as  one  green 
square. 


Figure  7.  Sensor  placement  by  (a)  Strategy  1  and  (b)  Strategy  2  when  the  obstacle  (black 
square)  coincides  with  the  high-value  area  (green  square).  Strategy  1  cannot  handle  the  case 
unless  the  high-value  area  is  wider  than  the  obstacle. 
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6  Safe-Mode  Concept 

The  safe-mode  concept,  presented  in  this  section,  is  an  application  of  the 
voting  logic  inference  (Klein  2007)  to  the  probabilistic  sensor  framework. 
Although  a  more  natural  application  of  this  inference  would  imply  a  sensor 
performance  measure  other  than  the  probability  of  detection  (e.g.,  a  de¬ 
gree  of  matching  of  a  received  signal  to  the  ideal  target  signature  or  num¬ 
ber  of  matching  signatures  for  multisignature  targets),  it  still  can  be  ap¬ 
plied  to  this  one. 

Sensor  performance  depends  on  operational  conditions  (e.g.,  combina¬ 
tions  of  terrain,  weather,  and  atmospheric  conditions)  at  a  specific  loca¬ 
tion.  For  example,  performance  of  the  acoustic  sensor  depends  substan¬ 
tially  on  wind  and  time  of  day  (Pettit  and  Wilson  2007).  In  other  words, 
the  radius  of  visibility  is  not  static,  it  depends  on  environmental  condi¬ 
tions.  The  most  conservative  way  to  account  for  this  performance  variabil¬ 
ity  is  to  always  use  the  sensor's  worst  anticipated  characteristics  in  a 
placement  strategy.  However,  operational  conditions  that  are  worst  for 
one  type  of  sensor  may  be  best  for  another,  so  that  a  situation  when  all  of 
the  sensors  have  their  worst  performance  simultaneously  may  never  occur. 
For  example,  visible-spectrum  cameras  are  unlikely  to  see  a  distant  object 
at  night  while  an  acoustic  sensor  has  the  highest  performance  at  this  time. 
The  safe-mode  concept  can  help  to  overcome  this  issue. 

Suppose  that  the  required  probability  of  detection  of  a  sensor  network  is 
equal  to  or  greater  than  0.9  (at  some  fixed  probability  of  false  alarm).  In 
the  case  of  two  independent  sensors,  there  are  three  modes  when  this  re¬ 
quirement  is  satisfied:  either  one  or  the  other  sensor  has  a  probability  of 
detection  over  0.9 ,  or  both  of  them  have  a  probability  of  detection  in  the 
range  0.7  -  0.9 .  Each  of  these  modes  is  safe  because  it  matches  the  total 
probability  of  detection  required.  Thus,  a  safe  mode  can  be  specified  by  a 
requirement  that  the  individual  sensor's  probability  of  detection  belong  to 
certain  bands  in  the  probability  space.  For  example,  Table  1  demonstrates 
these  bands  for  three  independent  sensors,  and  the  definition  of  the  safe 
modes  for  this  case  is  given  in  Table  2. 
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Table  1.  Band  definitions  for  three  independent  sensors. 


Band 

Pd 

Pfa 

1 

Uncertain  detection 

0-0.55 

10-6 

2 

Low  confidence 

0.55-0.7 

10-6 

3 

Medium  confidence 

0.7-0. 9 

10-6 

4 

High  confidence 

0.9-1 

10-6 

Table  2.  Definition  of  safe  modes  for  the  case  shown  in  Table  1. 


Safe 

Sensor 

Mode 

1 

2 

3 

1 

Band  2 

Band  2 

Band  2 

2 

Band  3 

Band  3 

3 

Band  3 

Band  3 

4 

Band  3 

Band  3 

5 

Band  4 

6 

Band  4 

7 

Band  4 

Each  of  the  safe  modes  provides 
a  probability  of  detection  greater 
than  0.9. 


Blank  spaces  signify  any  probability  of 
detection  value. 


As  operational  conditions  vary,  the  sensors'  probabilities  of  detection  also 
vary.  Suppose  there  are  H  combinations  of  weather/terrain  conditions  for 
which  Pd  of  each  sensor  has  been  calculated.  (In  practice,  this  could  be  a 

tedious  task).  This  H  is  assumed  to  be  large  enough  to  give  a  robust  statis¬ 
tic  about  operational  conditions.  Having  a  distribution  of  detection  prob¬ 
abilities  Pd ,  one  can  calculate  the  probability  that  a  sensor's  Pd  belongs  to 

a  certain  band.  Figure  8  shows  an  example  of  such  calculations  for  three 
independent  sensors.  Using  this  information,  the  main  question  could  be 
addressed:  what  is  the  probability  that  at  least  one  of  the  safe  modes  oc¬ 
curs  (with  respect  to  possible  operational  conditions)? 

To  answer  this  question,  one  should  calculate  the  probability  of  any  of  the 
safe  modes  that  are  not  mutually  independent  or  exclusive  (the  bands  for 
the  same  sensor  are  exclusive  because  they  are  chosen  with  no  overlap¬ 
ping,  but  the  modes  are  not).  The  answer  also  depends  on  the  location, 
meaning  that  the  safe-mode  probability  is  different  for  different  spatial 
points.  This  makes  it  possible  to  use  the  safe-mode  probability  in  a  place¬ 
ment  strategy  instead  of  the  probability  of  detection.  The  optimal 
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1  2  3 

Sensor 


Figure  8.  Probability  that  the  probability  of  detection  of  a  single  sensor  falls  into  a  certain 
band.  The  statistics  is  obtained  from  100  random  numerical  realizations  of  signal  and  noise 

parameters. 


sensor  locations,  then,  should  maximize  the  safe-mode  probability,  which 
means  a  high  robustness  relative  to  arbitrary  weather/terrain  conditions. 
Figure  9  shows  an  example  of  how  an  average  safe-mode  probability  de¬ 
pends  on  the  number  of  sensors  in  a  network.  The  red  line  indicates  the 
desired  level  of  the  safe-mode  probability  (equal  to  0.95  in  this  example). 
Note  that  if  the  safe-mode  probability  equals  one,  then,  under  any  opera¬ 
tional  conditions,  at  least  one  safe  mode  will  occur. 

The  application  of  this  concept  to  a  real  network  of  sensors  could  be  lim¬ 
ited  by  computational  time.  The  main  problem  is  an  exponential  growth  of 
safe  modes  as  the  number  of  sensors  increases.  Figure  10  demonstrates 
the  assessment  of  the  computational  time  (on  an  Intel  Xeon  2.66  GHz  PC 
with  4  Gb  of  RAM).  The  plot  consists  of  the  actual  elapsed 
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Figure  9.  Average  probability  of  safe  mode  versus  number  of  sensors  in  a  network. 


time  (the  blue  curve  with  markers)  and  the  prediction  (the  red  curve  with 
no  markers)  derived  from  a  curve  fitting  model  (the  first  experimental 
point  was  not  included  in  the  model).  The  elapsed  computational  time  te  is 

proportional  to  exp[M15],  where  M  is  the  total  number  of  sensors  in  a 
network.  According  to  Figure  to,  one  would  need  approximately  eight  full 
days  to  process  eight  sensors.  A  code  optimization  may  slightly  reduce  the 
time  required,  but  the  fact  that  the  safe  modes  are  not  independent  pre¬ 
vents  effective  calculations  of  the  safe-mode  probability.  The  practical  ap¬ 
plications  of  this  concept  directly  are  limited  to  seven  to  eight  sensors.  One 
way  to  incorporate  this  concept  for  a  larger  number  of  sensors  is  to  place 
sensors  in  accordance  with  the  probability  of  detection,  as  described  in 
sections  3  and  4,  and  then  model  how  this  probability  changes  due  to  vari¬ 
ous  operational  conditions  (already  at  fixed  sensor  locations).  Note  that 
this  approach  can  assess  the  safe-mode  probability  for  a  given  sensor's 
configuration  but  cannot  derive  the  most  robust  locations. 
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Total  number  of  sensors 


Figure  10.  Actual  (blue  curve  with  markers)  and  predicted  (red  curve)  computational  times 
required  for  safe-mode  solutions  versus  the  total  number  of  sensors. 
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7  Summary  and  Conclusions 

This  technical  report  presents  several  results  regarding  an  optimal  sensor 
placement.  Sensor  performance  is  evaluated  in  terms  of  the  probability  of 
detection  (at  some  given  probability  of  false  alarm).  A  physical  model  of 
signal  propagation  from  a  source  to  a  sensor  is  developed,  allowing  one  to 
obtain  a  realistic  spatial  dependence  of  the  probability  of  detection.  It  is 
unlikely  that  this  dependence,  shown  in  Figure  2,  can  be  described  by  the 
exponential  or  Gaussian  functions. 

Two  strategies  for  optimal  sensor  placement  were  presented.  Strategy  1  is 
heuristic  and  fast,  and  Strategy  2  is  a  solution  of  the  strict  binary  linear 
programming  problem.  Both  strategies  incorporate  the  coverage  prefer¬ 
ences  that  are  formulated  as  the  desired  probability  of  detection  and  can 
be  nonuniform  (allowing  one  to  assign  the  higher  probability  of  detection 
for  high-value  objects).  Strategy  1  places  sensors  one-by-one  in  the  loca¬ 
tion  where  they  are  most  needed  (the  worst  covered  and  the  farthest  away 
from  other  sensors)  while  Strategy  2,  strictly  speaking,  does  not  allow  the 
sequential  placement  and  yields  the  global  optimum  of  total  costs  associ¬ 
ated  with  sensors.  These  costs  can  include  not  only  the  number  of  sensors 
but  also  other  disincentives,  such  as  the  actual  price  of  sensors,  their  vul¬ 
nerability,  or  their  life-cycle  costs.  If  the  spatial  resolution  is  fine,  a  strict 
solution  of  Strategy  2  becomes  impractical  because  of  the  computational 
time  required.  For  this  case,  a  fast  algorithm  for  an  approximate  solution 
was  developed. 

To  compare  these  two  strategies,  a  numerical  experiment  was  conducted. 
In  this  experiment,  several  cases  were  studied,  such  as  the  coverage  of  a 
uniform  area,  an  area  with  high-value  objects,  and  an  area  with  obstacles. 
The  most  complicated  case,  yet  of  great  practical  importance,  was  also 
studied,  where  the  obstacles  coincide  with  high-value  objects.  The  numeri¬ 
cal  experiment  revealed  that  the  approximate  solution  of  the  strict  optimi¬ 
zation  problem  (Strategy  2)  outperforms  a  strict  solution  of  the  approxi¬ 
mate  optimization  problem  (Strategy  1). 

The  safe-mode  concept,  described  in  this  report,  views  the  sensor  place¬ 
ment  problem  from  the  standpoint  of  robustness  relative  to  different  op- 
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erational  conditions.  Operational  conditions  are  defined  by  combinations 
of  terrain,  weather,  and  daytime  conditions  in  the  place  where  a  sensor  is 
located.  Potentially,  this  concept  allows  one  to  find  the  most  robust  sensor 
placement.  However,  direct  application  of  this  concept  depends  signifi¬ 
cantly  on  the  total  number  of  sensors  and  is,  in  practice,  limited  to  several 
sensors.  One  way  to  bypass  this  limitation  is  to  use  worst-case  characteris¬ 
tics  of  each  sensor  in  the  placement  strategies  accounting  for  the  probabil¬ 
ity  of  detection.  This  will  provide  a  high  value  of  the  safe-mode  probability, 
although  may  be  excessive  because  the  situation  when  all  of  the  sensors 
have  their  worst  characteristics  may  never  occur.  Another  way  to  adopt  the 
safe-mode  concept  to  a  greater  number  of  sensors  is  to  employ  distributed 
calculations  (under  condition  that  an  algorithm  for  safe  modes  can  be  par¬ 
allelized).  This  possibility  was  not  studied  in  this  technical  report. 


ERDC/CRREL  TR-08-24 


36 


References 


Balanis,  C.  A.  1989.  Advanced  engineering  electromagnetics.  New  York:  J.  Wiley  &  Sons. 

Burdic,  W.  S.  1984.  Underwater  acoustic  system  analysis.  Englewood  Cliffs,  NJ: 
Prentice-Hall. 

Dhillon,  S.  S.,  and  K.  Chakrabarty.  2003.  Sensor  placement  for  effective  coverage  and 
surveillance  in  distributed  sensor  networks.  Proc  IEEE  Wireless 
Communications  and  Networking  Conference  3:1609-1614. 

Eidenbenz,  S.  2002.  Approximation  algorithms  for  terrain  guarding.  Inf  Process  Lett 
82:99-105. 

Elnagar,  A.,  and  L.  Lulu.  2005.  An  art  gallery-based  approach  to  autonomous  robot 
motion  planning  in  global  environments.  Proc  IEEE/RSJ International 
Conference  on  Intelligent  Robots  and  Systems  2079-2084. 

Erdem,  U.  M.,  and  S.  Sclaroff.  2004.  Optimal  placement  of  cameras  in  floor  plans  to 

satisfy  task  requirements  and  cost  constraints.  In  Proceedings,  5th  Workshop  on 
Omnidirectional  Vision,  Camera  Networks  and  Non-Classical  Cameras,  May  16, 
2004.  Prague,  Czech  Republic. 

Gonzalez-Banos,  H.,  A.  Efrat,  J.  C.  Latombe,  E.  Mao,  and  T.  M.  Murali.  1999.  Planning 
robot  motion  strategies  for  efficient  model  construction.  In  Proceedings,  gth 
International  Symposium  on  Robotics  Research.  Salt  Lake  City,  UT. 

Ishimaru,  A.  1978.  Wave  propagation  and  scattering  in  random  media.  New  York: 
Academic  Press.  (Reissued  1997,  New  York:  IEEE/OUP). 

Klein,  L.  A.  2007.  Sensor  and  data  fusion:  a  tool  for  information  assessment  and 
decision  making.  Washington,  DC:  SPIE  Press. 

McKenna,  J.  1987.  Worst-case  optimal  hidden-surface  removal.  ACM  Trans  Graphics 
6:19-28. 

Mittal,  A.,  and  L.  S.  Davis.  2008.  A  general  method  for  sensor  planning  in  multi-sensor 
systems:  extension  to  random  occlusion.  IntJ  Computer  Vision  76(i):3i-52. 

O’Rourke,  J.  1987.  Art  gallery  theorems  and  algorithms.  New  York:  Oxford  University 
Press. 

Ostashev,  V.  E.  1997.  Acoustics  in  moving  inhomogeneous  media.  London,  UK:  E&FN 
Spon. 

Pettit,  C.  L.,  and  D.  K.  Wilson.  2007.  Proper  orthogonal  decomposition  and  cluster 
weighted  modeling  for  sensitivity  analysis  of  sound  propagation  in  the 
atmospheric  surface  layer.  J AcoustSocAm  122(3):1374-1390. 


ERDC/CRREL  TR-08-24 


37 


Shermer,  T.  C.  1992.  Recent  results  in  art  galleries.  ProcIEEE  8o(9):i384-i399. 

Sierksma,  G.  2002.  Linear  and  integer  programming:  theory  and  practice,  2nd  Ed.  New 
York;  Basel,  Switzerland:  Marcel  Dekker. 

Tatarskii,  V.  1. 1971.  The  effects  of  the  turbulent  atmosphere  on  wave  propagation. 
Jerusalem,  Israel:  Keter  Publishing  House. 

Wilson,  D.  K.,  C.  L.  Pettit,  M.  S.  Lewis,  S.  Mackay,  and  P.  Seman.  2008.  Probabilistic 
framework  for  characterizing  uncertainty  in  the  performance  of  networked 
battlefield  sensors.  Proc  SPIE  698104:1-12. 

Wixson,  L.  1994.  Viewpoint  selection  for  visual  search.  Proc  IEEE j Conference  on 
Computer  Vision  and  Pattern  Recognition  800-805. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining 
the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for 
reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington, 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not 
display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

December  2008  Technical  Report 


4.  TITLE  AND  SUBTITLE 

Optimal  Sensor  Placement  with  Terrain-Based  Constraints  and 
Signal  Propagation  Effects 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Sergey  N.  Vecherin,  D.  Keith  Wilson,  and  Chris  L.  Pettit 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Engineer  Research  and  Development  Center 
Cold  Regions  Research  and  Engineering  Laboratory 
72  Lyme  Road 
Hanover,  NH  03755-1290 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


ERDC/CRREL  TR-08-24 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

The  optimal  sensor  placement  problem,  as  considered  here,  is  to  select  the  types  and  locations  of 
sensors  providing  coverage  at  high-value  terrain  locations  while  minimizing  a  specified  cost  func¬ 
tion.  The  cost  function  can  reflect  various  disincentives,  such  as  the  actual  cost  of  the  sensors, 
the  total  number  of  sensors,  and  the  probability  that  the  sensor  will  be  found  and  disabled  by  hos¬ 
tile  actors.  The  probability  of  detection  (at  a  certain  probability  of  false  alarm)  is  assumed  to 
depend  on  terrain  conditions  and  obstructions,  and  may  be  arbitrarily  complex.  Two  strategies  are 
described  for  finding  the  minimal  number  of  sensors,  and  their  locations  that  will  satisfy  given  cov¬ 
erage  preferences.  The  first  is  heuristic  in  nature  and  based  on  placing  sensors  one-by-one  where  the 
probability  of  detection  is  minimal.  This  strategy  offers  a  rapid,  but  suboptimal  solution.  The  sec¬ 
ond  strategy  is  based  on  solution  of  the  binary  linear  programming  problem.  For  the  case  of  fine  spa¬ 
tial  resolution  that  leads  to  large  matrix  dimensions,  a  fast  algorithm  for  approximate  solution  of 
this  problem  is  developed.  The  key  features  of  this  study  are:  1)  the  probabilistic  framework  of  sen¬ 
sor  performance,  2)  incorporation  of  the  coverage  preferences  in  the  placement  strategy,  3)  realistic 
modeling  and  incorporation  of  the  sensors'  probability  of  detection,  4)  multimodal  sensor  support,  5) 
a  strict  formulation  of  the  optimal  coverage  problem,  6)  development  of  a  fast  algorithm  for  approxi¬ 
mate  solution  of  the  binary  linear  programming  problem,  and  7)  introduction  of  a  safe-mode  concept. 


15.  SUBJECT  TERMS 

binary  linear  programming,  multimodal  sensor  network,  optimal  sensor  placement,  probabilis¬ 
tic  framework,  probability  of  detection,  safe-mode  concept,  sensor  location  optimization, 
sensor  performance,  signature  propagation  effects,  weather  effects  on  sensors 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE 

OF  ABSTRACT 

OF  PAGES 

PERSON 

a.  REPORT 

U 

b.  ABSTRACT 

u 

c.  THIS  PAGE 

U 

U 

45 

19b.  TELEPHONE  NUMBER  (include 
area  code) 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


